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ABSTRACT 

Conventional dynamic model of gamma-ray burst remnants is found to be 
incorrect for adiabatic blastwaves during the non-relativistic phase. A new 
model is derived, which is shown to be correct for both radiative and adiabatic 
blastwaves during both ultra-relativistic and non-relativistic phase. Our model 
also takes the evolution of the radiative efficiency into account. The importance 
of the transition from the ultra-relativistic phase to the non-relativistic phase is 
stressed. 

PACS: 95.30.Lz, 98.70.Rz, 97.60. Jd 

The origin of gamma-ray bursts (GRBs) has remained unknown for over 30 years. 1 ' 2 A 
major breakthrough appeared in early 1997, when the Italian-Dutch BeppoSAX satellite 
observed X-ray afterglows from GRB 970228 for the first time. 3 Since then, X-ray afterglows 
have been observed from about 15 GRBs, of which ten events were detected optically and 
five bursts were also detected in radio wavelengths. The cosmological origin of at least some 
GRBs is firmly established. The so called fireball model 4 ' 5 is strongly favoured, which is 
found successful at explaining the major features of the low energy light curves. 6 ~ 9 

In the fireball model, low energy afterglows are generated by ultra-relativistic fireballs, 
which first give birth to GRBs through internal or external shock waves and then decelerate 
continuously due to collisions with the interstellar medium (ISM). The dynamics of the 
expansion has been investigated extensively 6-9 Both analytic solutions and numerical 
approaches are available. It is a general conception that current models describe the gross 
features of the process very well and further improvements are possible only by considering 
some details. However, we find that three serious problems are associated with the popular 
model. 
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First, it is usually assumed that the expansion is ultra-relativistic. Then for an 
adiabatic fireball, the evolution of the bulk Lorentz factor is derived to be: 

7 » (200 - 400) El( 8 ni 1/8 r 3/8 , (1) 

where £51 = -Eo/(10 51 erg) with E the initial fireball energy, n\ = n/(lcm -3 ) with n the ISM 
number density, and t is observer's time in unit of s. 6-9 The radius of the blastwave scales 
as R oc t 1 / 4 . Based on Eq.(l), flux density at frequency v then declines as S u oc t 3 ^ 1_p ^ 4 , 
where p is the index characterizing the power-law distribution of the shocked ISM electrons, 
dn' e /d^ c oc 7~ p . These expressions are valid only when 7>1. 

However, optical afterglows from GRB 970228 and GRB 970508 were detected for 
as long as 190 and 260 d respectively, while in Eq.(l), even t = 30 d will lead to 7 ~ 1. 
It is clear that the overall evolution of the postburst fireball can not be regarded as 
a simple one-phase process, we should pay special attention to the transition from the 
ultra-relativistic phase to the non-relativistic phase. 10 This is unfortunately ignored in the 
literature. 

Second, the expansion of the fireball might be either adiabatic or highly radiative. 
Extensive attempts have been made to find a common model applicable for both cases. 11-13 
As a result, a differential equation has been proposed by various authors, 12 ' 13 

dm M ' K J 

where m is the rest mass of the swept-up ISM, M is the total mass in the co-moving 
frame, including internal energy U. Since thermal energy produced during the collisions is 
dE = c 2 (7 — l)dm, usually we assume: 13 

dM = ' ~ € > dE + dm= [(1 - e) 7 + e]dm, (3) 

where e is defined as the fraction of the shock generated thermal energy (in the co-moving 
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frame) that is radiated. It is putative that Eq.(2) is correct in both ultra- relativistic and 
non-relativistic phase, for both radiative and adiabatic fireballs. 

In the highly radiative case, e = 1, dM = dm, Eq.(2) reduces to, 

c?7 7 2 — 1 



dm M C j + m ' 

where M e j is the mass ejected from the GRB central engine. Then an analytic solution is 

available, 11 ' 13 which satisfies 7 ex R~ 3 when 7>1 and v oc R~ 3 when 7 ~ 1, where v is 

the bulk velocity of the material. These scaling laws indicate that Eq.(2) is really correct 

for highly radiative fireballs. In the adiabatic case, e = 0, dM = / ~fdm, Eq.(2) also has an 

analytic solution: 12 

M = [M 2 + 2 7o M cj m + m 2 ] 1/2 , (5) 

m + 7o M ej 
7 M ' W 

where 70 is the initial value of 7. During the ultra-relativistic phase, Eqs.(5) and (6) do 

produce the familiar power-law 7 oc i?~ 3//2 , which is often quoted for an adiabatic blastwave 

decelerating in a uniform medium. In the non-relativistic limit (7 ~ 1, m 3> 7oM e j), 

Chiang and Dermer have derived 7^1 + 7 M C j/m, 12 so that they believe it also agrees 

with the Sedov solution, i.e., v ex i? -3 / 2 . 14 However we find that their approximation is 

not accurate, 15 because they have omitted some first-order infinitesimals of 7oM C j/m. The 

correct approximation could be obtained only by retaining all the first and second order 

infinitesimals, which in fact gives: 7^1 + (7 M C j/m) 2 /2, then we have v oc i?~ 3 . 15 This is 

not consistent with the Sedov solution! 

The problem is serious: (i) It means that the reliability of Eq.(2) is questionable, 
although it does correctly reproduce the major features for radiative fireballs and even 
for adiabatic fireballs in the ultra-relativistic limit, (ii) In the non-relativistic phase of 
the expansion, the fireball is more likely to be adiabatic rather than highly radiative. 
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However, it is just in this condition that the conventional model fails. So any calculation 
made according to Eq.(2) will lead to serious deviations in the light curves during the 
non-relativistic phase. 

Third, for simplicity, it is usually assumed that e is a constant during the expansion. 
But in realistic case this is not true. The fireball is expected to be highly radiative (e = 1) 
at first, due to significant synchrotron radiation. In only one or two days, it will evolve to 
an adiabatic one (e = 0) gradually. So e should evolve with time. 16 

Below, we will construct a new model that is no longer subject to the aforementioned 
problems. 

In the fixed frame, since the total kinetic energy of the fireball is 
Er = (7 _ l)(-^ej +'m)c 2 + (1 — e)7?7, 17 and the radiated thermal energy is €7(7 — l)dmc 2 , n 
we have: 

d[(7 - l)(M cj + m)c 2 + (1 - e) 7 £7] = -e 7 ( 7 - l)c 2 dm. (7) 

For the item U, it is usually assumed: dll = c 2 (7 — l)dm. 17 Eq.(2) has been derived 

just in this way. However, the jump conditions 11 at the forward shock imply that 

U — (7 — l)mc 2 , so we suggest that the correct expression for dU should be: 

dU = ^[(7 — l)mc 2 } = (7 — l)dm c 2 + mc 2 d r ). Here we simply use U — (7 — l)mc 2 and 

substitute it into Eq.(7), then it is easy to get: 15 

^7 = 7 2 ~1 (g) 

dm M ej + em + 2(1 - e)^m 

In the highly radiative case (e = 1), Eq.(8) reduces to Eq.(4) exactly. While in the 
adiabatic case (e = 0), Eq.(8) reduces to: 

*L = l2 ~ l . (9) 

dm M e j + 27m 

The analytic solution is: 

(7 - l)M ej c 2 + ( 7 2 - l)mc 2 = E . (10) 
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Then in the ultra-relativistic limit, we get the familiar relation of 7 oc i?~ 3 / 2 , and in the 
non-relativistic limit, we get v oc R~ 3 ^ 2 as required by the Sedov solution. From these 
analysises, we believe that Eq.(8) is really correct for both radiative and adiabatic fireballs, 
and in both ultra-relativistic and non-relativistic phase. 

In realistic fireballs, e is a variable dependent on the ratio of synchrotron-radiation- 
induced to expansion-induced loss rate of energy 16 As usual, we assume that in the 
co-moving frame the magnetic field energy density is a fraction £| of the energy 
density e', B' 2 /(87r) = £|e', and that the electron carries a fraction £ e of the energy, 
7c,min = £c(7 — l)m p /m e + 1, where m p and m c are proton and electron masses, respectively. 
The co- moving frame expansion time is t' cx = R/^c), and the synchrotron cooling time is 
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t' syn = 67rm c c/(o- T £> /2 7 miniC ), where a T is the Thompson cross section. Then we have 

t'- 1 
e = £ c / yn - . (11) 

syn ' "ex 

We have evaluated Eqs.(8) and (11) numerically, bearing in mind that: 18 

dm = 47iR 2 nm p dR, (12) 



dR = w 7 ( 7 + yV - l)dt. (13) 

We take E = 10 52 erg, n = 1 cm" 3 , M cj = 2 x 10" 5 M . Figs.(l) and (2) illustrate the 
evolution of 7 and R, where full, dotted, and dashed lines correspond to constant e values 
of 0, 0.5, and 1 respectively. Dash-dotted lines are plot by allowing e to vary according to 
Eq.(ll). It is clearly shown that our new model overcomes the shortcomings of Eq.(2). For 
example, for highly radiative expansion, the dashed lines in these figures approximately 
satisfy 7 oc t~ 3 ^ 7 , R oc t 1 / 7 , 7 oc R~ 3 when 7 ^> 1, and v oc t _3//4 , R oc t 1 / 4 , v oc R~ 3 when 
7 ~ 1. While for adiabatic expansion, the full lines satisfy 7 oc t~~ 3 ' 8 , R oc t 1 ' 4 , 7 oc R~ 3 ' 2 
when 7 ^> 1, and satisfy v oc £ -3 / 5 , R oc t 2//5 , v oc R~ 3 / 2 when 7 ~ 1. 
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In order to compare with observations, we have also calculated the synchrotron 
radiation from the shocked ISM. Fig. (3) illustrates R band afterglows from the realistic 
fireball. We see that after entering the non-relativistic phase, the light curve becomes 
steeper only slightly, consistent with the prediction made by Wijers et al.? In contrast, 
Eq.(2) generally leads to a much sharper decline. 18 In our new model optical afterglows 
from GRB 970228 are generally well fitted. 

To conclude, current researches are mainly concentrated on the ultra-relativistic phase 
of the expansion of GRB remnants. The popular dynamic model is in fact incorrect for 
adiabatic fireballs during the non-relativistic phase. This is completely unnoticed in the 
literature. We have revised the model. Our new model has been shown to be correct in 
both ultra-relativistic and non-relativistic phase. The revision is of great importance, taking 
account of the following facts: (i) Optical afterglows lasting for more than 100 - 200 d have 
been observed from some GRBs, the advent of the non-relativistic phase seems inevitable, 
(ii) Beaming effects also lead to a steepening in the optical light curve, non-relativistic 
effects should be considered carefully to tell whether GRB ejecta are beamed or not, which 
is crucial in understanding the GRB origin, (iii) HI supershells might be highly evolved 
GRB remnants, 19 ' 20 to address this question in detail, we should deal with non-relativistic 
blastwaves. Additionally we suggest that at very late stages, GRB remnants might become 
highly radiative again, just in the same way that supernova remnants do. 14 This might 
occur when the bulk velocity is just several tens kilometers per second. 
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Figure Caption 

Fig. 1. — Evolution of the bulk Lorentz factor 7. We take E$ = 10 52 erg, n = 1 cm~ 3 , 
M e j = 2 x 10~ 5 M Q . The full, dotted, and dashed lines correspond to e = (adiabatic), 0.5 
(partially radiative), and 1 (highly radiative) respectively. The dash-dotted line is plotted 
by allowing e to evolve with time. 

Fig. 2. — Evolution of the radius. Parameters and line styles are the same as in Fig.l. 

Fig. 3. — R band afterglows from a realistic fireball. We take p = 2.1, £ c = 1.0, £| = 0.01. 
The distance D is 18 Gpc. Other parameters are the same as in Fig.l. For the origin of the 
observational data points, please see Ref. [10]. 



This manuscript was prepared with the AAS IATpX. macros v4.0. 




J — I — I — lo 






CV2 



CV2 

I 



I 



CO 

I 



(I -I) S T 



\l 1 1 1 1 1 1 1 1 1 1 




- \\ 


- 


\\ 




- \\ 


— 


\\ 




\\ 


_ 


\n \ \ 




■ \\ \ 


- 


\\ \ \ 




\\ \ \ 




\\ \ \ 




\\ \ \ 


- 


\\ \ \ 






_ 


\\\ \ 




%\ 


- 


\\\ \ 




\\\ \ 




\ v *. \ 




\ '*• \ 


- 




- 


\ w 




— \ " \ 

\ ,J A 




V* 


- 


Y>\ 




\* 


- 


\ 


- 


■ i i 1 i i i 1 i i A 


- 



_ in 



_ o 






- in 



CV2 
CV2 



O 
C\2 



CO 



CO 



(uio) h §i 



i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — tt| — i — r 

/' 



/ 



- Z> 



/ 



i.nt 



/* 



/' I — 4. 



/ 



/ 



/ M 1 

h^n 1 



/ 



\-M — 1 



■U 



- CO 






- 10 



00 

CV2 

O 

z> 

05 

PQ 



- ^ 



CV2 



J Ll I I I I I I I I I I I I I I I I I I L 

iH O *h CV2 

I I 



(A[H) * S §1 



